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Abstract 

Pb-2.2 wt% Sb alloy was directionally solidified in 1,2,3 and 7 mm diameter crucibles. 
Pb-Sb alloy presents a solutally unstable case. Under plane-front conditions, the resulting 
macrosegregation along the solidified length indicates that convection persists even in the 1 mm 
diameter crucible. Al-2 wt% Cu alloy was directionally solidified because this alloy was 
expected to be stable with respect to convection. Nevertheless, the resulting macrosegregation 
pattern and the microstructure in solidified examples indicated the presence of convection. 
Simulations performed for both alloys show that convection persists for crucibles as small as 0.6 
mm of diameter. For the solutally stable alloy, Al-2 wt% Cu, the simulations indicate that the 
convection arises from a lateral temperature gradient. 
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1. Introduction 


It is well established that convection influences the planar to cellular transition [1-4] and 
primary cellular/dendritic spacings [5-7]. During directional solidification the thermal profile in 
the melt provides stability against convection because the melt density decreases with increasing 
temperature. However, the solutal profile in the melt is stabilizing only for those alloys where 
the increased solute content results in higher melt density ( e.g ., hypoeutectic Al-Cu alloys). Even 
for these alloys the radial thermal gradient at the triple junction of the solid-liquid interface and 
the ampoule wall causes mild convection that is localized in the vicinity of the liquid-solid 
interface [8-9]. The solutal profile in the melt is destabilizing for those alloys in which the melt 
density decreases with increased solute content (e.g., hypoeutectic Pb-Sb alloys). It causes 
extensive convection in the overlying and interdendritic melt and results in macrosegregation 
along the solidified length for the planar, cellular or dendritic morphologies [10-11]. 

Since the theoretical dendrite growth and morphology models consider only the diffusive 
thermal and solutal transports, attempts have been made to eliminate convection in order to allow 
for a quantitative comparison of the theoretical predictions with the experimental observations 
[9,12]. Effecting directional solidification without convection could also provide data that could 
be used to deduce diffusion coefficients in alloy melts. In a recent study it has been shown that 
reducing the ampoule diameter to about 1 mm almost completely eliminates convection during 
directional solidification of Al-4 wt% Cu alloy [8,9]. The purpose of this research was to 
examine if a similar reduction in the ampoule diameter can reduce the more intense, longer 
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range, convection that occurs during directional solidification of Pb-2.2 wt% Sb alloy. The 
experiments are supported by simulations of the convection, using a finite element model that 
has been used for the solidification of alloys [13-15]. 


2. Experimental procedures 

Since the directional solidification furnace is described elsewhere [10-1 1], only a brief 
description of its operation is provided here. The furnace assembly is withdrawn and the sample 
is held stationary during directional solidification. In addition, the whole assembly sits on top of 
a vibration isolation platform. Samples are kept in quartz ampoules and are directionally 
solidified in a flowing argon atmosphere. The furnace is initially translated downwards to remelt 
the cast alloy feed stock and about 1 cm of a 5 cm long seed crystal (length parallel to [100]). 

The furnace is then translated upwards to achieve the directional solidification. After solidifying 
for about 10 cm, the rest of the melt column is quenched by withdrawing the furnace rapidly and 
blasting the ampoule surface with helium gas that is cooled in a tank of liquid nitrogen. 

Cast Pb-2.2 wt% Sb feedstock-rods were placed into quartz ampoules (7 mm inner 
diameter) on top of the single crystal seeds. An ampoule also contained a bundle of 1 , 2, and 3 
mm inner diameter quartz capillaries that were connected to a vacuum. The furnace was first 
lowered in order to melt the charge and about 1 cm portion of the seed. A 20-22 cm long melt- 
column was then sucked into the capillaries by inserting the capillary bundle into the melt and 
activating the vacuum. The capillaries were further lowered till their bottoms were near the 
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liquid-solid interface, and then the entire assembly was directionally solidified at 10 pm s' 1 for a 
distance of about 5 cm. This method of casting provided a directionally solidified portion at the 
bottom of the capillaries that was parallel to [100] and ensured that during subsequent remelting 
and directional solidification the samples grew along [100]. 

The pre-cast 1 , 2 and 3 mm ID capillary samples were bundled together and kept inside 
another 7 mm ID quartz ampoule. The capillary bundle was surrounded by gallium metal kept in 
the 7 mm ED ampoule in order to ensure that all of the samples solidified under identical thermal 
conditions. After remelting about 14 cm long portion (leaving about 3 cm long [100] oriented 
solid at the bottom), the capillary samples were directionally solidified at 0.4 pm s’ 1 at a thermal 
gradient of 82 K cm' 1 . Two 7 mm diameter samples were also directionally solidified at 120 
K cm' 1 and a growth speed of 0.4 pm s' 1 in order to examine the influence of transverse magnetic 
field (0.45 T) on convection. One sample was grown in the transverse magnetic field, and the 
other without the field. 

Samples with dendritic morphologies were obtained by directional solidification at 3 and 
30 pm s' 1 with a thermal gradient of 82 K cm' 1 . The growth speed of 3 pm s' 1 corresponds to the 
dendrites in the vicinity of the cell to dendrite transition, whereas at 30 pm s' 1 the morphology is 
within the fully dendritic regime. 

The directionally solidified samples were sectioned along their lengths in order to obtain 
slices, about 2 mm thick, for chemical analysis by atomic absorption spectroscopy. A specimen 
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about 2 cm long in the vicinity of the quenched liquid-solid interface was used for optical 
metallography. 

In addition to directional solidification (DS) of the Pb-Sb alloy, a DS experiment on 
Al-2 wt% Cu alloy was also performed. The Al-Cu alloy was solidified in an alumina tube with 
an internal diameter of 6 mm. This alloy was solidified in a thermal gradient of 90 K cm' 1 at 0.5 
pm s' 1 . It was also solidified in the transverse magnetic field (0.45 T) and with no magnetic 
field. 


3. Experimental results 

Directional solidification of Pb-2.2 wt% Sb results in thermally stable and solutally 
unstable density profile in the melt above the planar liquid-solid interface. This results in 
extensive convection of the liquid as indicated by the experimentally observed longitudinal 
macrosegregation and computer simulations. For this alloy the planar liquid-solid interface is 
essentially flat with some curvature next to the quartz ampoule (Fig. 1(a)). On the other hand, 
the solutal profile and the thermal profile in the melt are expected to be stable against natural 
convection during directional solidification of Al-2 wt% Cu alloy, because copper enrichment 
results in increased melt density. Ideally one should expect a diffusive solutal profile during 
plane front growth of Al-Cu alloy. However, this alloy grown at 90 K cm' 1 and 0.5 pm s' 1 also 
shows evidence of convection as indicated by the convex (towards the liquid) shape of the liquid- 
solid interface (Fig. 1(b)). The major contributing factor that makes the shape convex is the 
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solutal accumulation near the outer portion of the sample. This is driven by jf convection due to 
the radial thermal gradients that is localized in the melt in the immediate vicinity of the liquid- 
solid interface. This convection, however, produces less severe longitudinal macrosegregation 
compared to the Pb-2.2 wt% Sb alloy as shown below. 


3.1. Influence of Transverse Magnetic Field 

Fig. 2 shows the variation in solute content in the Pb-2.2 wt% Sb and Al-2 wt% Cu alloy sample 
grown in 7 mm ID quartz and 6 mm ID alumina ampoules, respectively, and with and without 
the application of magnetic field. It plots C s /C„, where C s is the solute content after directional 
solidification and C 0 is the initial solute content. The solutal profiles in this figure are clearly 
not those expected from diffusive mass-transport. The segregation patterns indicate extensive 
convection during directional solidification of the Pb-Sb alloy and even some convection in the 
Al-Cu alloy. For diffusive transport, one would expect a steady-state value of unity for C s IC„ 
after an initial transient, during which the solute content increases from k C,„ where k is the 
equilibrium partition ratio. The length of the transient would be expected to be D/R k [16], 
where D is the solutal diffusivity in the melt and R is the growth speed. For example, the value 
of D/Rk is 0.9 cm for the Pb-2.2 wt% Sb alloy grown at 0.4 |lm s' 1 . The extent of longitudinal 
macrosegregation in the Pb-Sb alloy is significantly higher than that in the Al-Cu alloy, and an 
application of a 0.4 T transverse magnetic field does not have any influence on the extent of the 
macrosegregation . 
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3.2. Effect of Reducing the Ampoule Diameter 


Fig. 3 compares the dimensionless concentration along the solidified length for the 1,2,3 
and 7 mm diameter Pb-2.2 wt% Sb samples that were directionally solidified with a planar 
liquid-solid interface. The extent of macrosegregation shows a systematic decrease with the 
decreasing ampoule diameter. This indicates that the intensity of convection decreases with 
decreasing ampoule diameter, but that even with a diameter as small as 1 mm convection 
persists. 


4. Computer simulations 

Details of the solidification/convection model and the computer program are given 
elsewhere [13-15] and not repeated here. In the momentum and continuity equations, the 
following assumptions are invoked: the liquid is Newtonian and the flow is laminar; the 
Boussinesq approximation is made in the buoyancy term of the momentum equation; the solid 
phase is stationary; and the densities of solid ( p s ) and liquid ( p t ) are different but constant. The 
mushy zone is treated as a porous medium in which the intergranular/ interdendritic liquid 
convects within a solidifying dendritic solid-network. The model comprises comprehensive 
equations for continuity, conservation of momentum, conservation of energy, and a conservation 
equation for each element in the alloy. The flow of the liquid in the mushy zone is dictated by 
the permeability, which is a function of fraction liquid and a characteristic length-scale of the 
dendritic structure ( e.g ., primary dendrite area spacing). As the fraction liquid approaches unity 
the permeability becomes infinity, and the transport conservation equations reduce to those of the 
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all-liquid region. Thus, there is a natural transition and coupling of the transport phenomena in 
the mushy zone and the all-liquid zone. The model also accounts for the partitioning of each 
alloy element between solid and liquid. The density of the liquid (be it in the mushy zone or in 
the all-liquid zone) depends on its constitution and temperature; hence the liquid density varies 
spatially and temporally and is subject to gravity-driven convection. The fraction of liquid 
adjusts according to the constitution and freezing temperature of the liquid; hence both 
solidification and remelting during solidification are possible. An extensive list of properties for 
each alloy is required for running the computer simulations. These are presented in Tables 1 and 
2 for the Pb-2.2 wt% Sb alloy and the Al-2 wt% Cu alloy, respectively. 

4.1. Simulations of Convection in Pb-2.2 wt% Sb Alloy 

The alloy of Pb-2.2 wt% Sb is directionally solidified by simulation in a rectangular mold 
of 200 mm in height and various widths; gravity acts downward. Initially, the mold is filled with 
an all-liquid alloy of the nominal composition in a stable vertical temperature gradient such that 
the bottom temperature is slightly above that of the liquidus temperature. A constant thermal 
gradient (8200 K m -1 ) is imposed at the top, and a constant cooling rate (0.0164 K s -1 ) is applied 
at the bottom, which results in a solidification rate of 2 Jim s -1 . Lateral walls are thermally 
insulated; liquid flows through the top boundary to compensate for the shrinkage. A no-slip 
condition is used for velocity at the bottom and along the two vertical boundaries, and a stress- 
free condition is used on the top boundary. Solute diffusion flux is set to zero on the boundaries. 

The experimental castings were solidified at 0.4 p,m s _1 under plane front conditions. The 
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simulations duplicate the thermal gradient used in the experiment, but the simulated results are 
for 2 pm s' 1 . As explained above, the computer program is designed for simulating dendritic 
solidification with a mushy-zone of finite thickness. Solidification rates of less than 2 pm s' 1 
were attempted, but the calculated results were unrealistic and discounted. Similar 
computational anomalies were encountered and discussed by Huang et al. [35]. 

The solidification and convection in the 2-mm wide domain of Pb-2.2 wt% Sb was 
simulated, and the results of the bottom half at 2500 s are shown in Figs. 4(a) and 4(b). Fig. 4(a) 
shows thermosolutal convection as velocity vectors and contour lines indicating the mushy zone. 
Fig. 4(b) shows the total concentration of Sb. During solidification, the Sb-enriched layer in the 
liquid above the dendrite tip makes the liquid unstable with respect to convection, and there are 
well-developed convection cells with a width of approximately half of the domain width. The 
maximum fluid velocity is 162 pm s' 1 . Also notice that the convection results in the intrusions 
into the mushy zone at the intersection of the top of the mushy zone and the vertical boundaries. 

To see the effect of width of domain on the convection, the solidification domains with 
widths of 1, 0.6 and 0.3 mm were solidified with the same thermal conditions as the domain with 
2-mm width. As shown in Figs. 5(a) and 5(b), there are still well-developed convection cells for 
the cases with 1 and 0.6 mm of widths. It is observed that the maximum fluid velocity is 
decreased to 102 and 56 pm s' 1 , as the width is reduced to 1 and 0.6 mm, respectively. Finally, 
no convection is observed at 0.3 mm of the domain width, and there are only small downward 
velocities due to shrinkage as shown in Fig. 5(c). These simulations certainly support the 
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evidence of thermosolutal convection during solidification in capillaries as small as 1 mm 
(Fig. 3). 

The thermosolutal convection in the Pb-2.2 wt% Sb alloy is driven by the Sb-enriched 
layer just above the top of the mushy zone [36]. Solidification at rates less than 2 pm s -1 would 
cause the enriched layer to have yet higher concentration of Sb. That thermosolutal convection 
exists when the solidification rate is 2 pm s -1 tells us that there is also thermosolutal convection 
at slower solidification rates, at least for widths of 0.6 mm or greater. 

4.2. Simulations of Convection in Al-2 wt% Cu Alloy 

The alloy of Al-2 wt% Cu is directionally solidified by simulation in a rectangular mold 
of 200 mm in height and 6 mm in width; gravity acts downward. A constant thermal gradient 
(9000 K nrf ') is imposed at the top, and a constant cooling rate (0.036 K s -1 ) is applied at the 
bottom that results in a solidification rate of 4 pm s _1 . The vertical surfaces were thermally 
insulated. Other boundary conditions are same as the case of Pb-Sb alloy described above. Since 
the Cu-enriched layer in the liquid above the dendrite tip is solutally stable, no convection was 
observed in the domain with a width of 6 mm. 

The Cu-concentration in the liquid at the dendrite tips is C, = 3.4 wt%. When the 
solidification rate becomes slower to reach planar solidification, the liquid concentration at the 
dendrite tips is as high as C 0 Ik (i.e., 1 1.6 wt% with C 0 = 2 wt% and = 0.173). Hence, by 
artificially decreasing the equilibrium partition ratio of Cu ( k ) from 0.173 to 0.0173 and the slope 
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of liquidus ( m L ) from 3.433 to 0.344 K (wt%)'\ the liquid concentration at the dendrite tips (C,) 
can be increased to 11.6 wt% in the simulations. However with thermally insulated vertical 
sides, there was still no thermosolutal convection simulated with the higher liquid concentrations 
at the dendrite tips. Hence, provided the vertical surfaces are thermally insulated, the Al-Cu 
alloy exhibits no convection. 

To observe the effect of lateral temperature gradient, a heat flux (5.4 x 10 4 J s' 1 m‘ 2 ) was 
imposed on the lateral wall instead of imposing the thermally insulated condition. The lateral 
heat flux was estimated from the experimental temperature profile and is approximately 5 % of 
the heat flux through the bottom surface of the domain. As shown in Figs. 6(a) and 6(b), the 
convection arises from lateral temperature gradient of the domain width of 6 mm, and the 
maximum fluid velocity is 133 pm s’ 1 . Similar calculations were done for smaller widths, and 
the convection persists for crucibles as small as 0.6 mm of diameter as a result of the lateral 
temperature gradient. 
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5. Conclusions 


Following conclusions can be drawn from this study in which Pb-2.2 wt% Sb and Al-2.2 
wt% Cu alloys were directionally solidified with planar liquid-solid interfaces. 

(1) Convection during directional solidification of Pb-2.2 wt% Sb alloy, with a planar liquid- 
solid interface, produces extensive longitudinal macrosegregation. Application of a 0.4 T 
transverse magnetic field has no effect on the extent of this convection or segregation. Reducing 
the ampoule diameter decreases the extent of convection. However, elimination of thermosolutal 
convection would require solidification in capillaries with diameters less than 600 |0.m, which 
would introduce significant experimental difficulties. 

(2) Convection driven by lateral thermal gradients occurred during the directional solidification 
of the Al-2 wt% Cu alloy. 

(3) Simulations performed for both alloys show that convection persists for crucibles as small as 
0.6 mm. For the solutally stable alloy, Al-2 wt% Cu, the simulations indicate that the convection 
arises from lateral temperature gradient. 
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Table 1 Thermodynamic and transport properties of Pb-2.2 wt% Sb alloy 
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* The code used for the simulations requires that a finite value of D s be used, so D s is made approximately 
four orders magnitude smaller than D t . 













four orders magnitude smaller than D t . 












Figure captions 


Fig. 1. Typical morphology of liquid-solid interface during plane front directional solidification 
of Pb-2.2 wt% Sb and Al-2.2 wt% Cu alloy. Notice the convex liquid-solid interface for the Al- 
Cu alloy as compared with nearly flat interface for the Pb-Sb alloy. The radial solute 
segregation has lead to the break-down of the planar front near the periphery of the Al-Cu 
sample, where cellular morphology has evolved. 

Fig. 2. Macrosegregation along directionally solidified length of Pb-2.2 wt% Sb alloy grown in 7 
mm ID quartz crucibles at 0.4 pm s' 1 and 120 K cm" 1 ; and of an Al-2 wt% Cu alloy grown in 6 
mm ID alumina crucibles at 90 K cm' 1 and 0.5 pm s' 1 . Each alloy was grown in the presence of a 
transverse magnetic field of 0.4 T (MF) and without it (NMF). 

Fig. 3. Influence of reducing the ampoule diameter on the longitudinal macrosegregation in 
directionally solidified Pb- 2.2 wt% Sb alloy samples (plane front). The 1, 2, and 3 mm 
samples were grown at 0.4 pm s' 1 , 82 K cm' 1 , while the 7 mm sample was grown at 0.4 pm s' 1 , 
120 K cm' 1 . 

Fig. 4. Simulated results of Pb-2.2 wt% Sb at 2500 s: (a) thermosolutal convection and contour 
lines indicating mushy zone; (b) total concentration of Sb, wt%. Dimensions of domain are 
in m. 

Fig. 5. Simulated results of Pb-2.2 wt% Sb at 2500 s: domain widths of (a) 1 mm; (b) 0.6 mm; 
and (c) 0.3 mm. Contour lines indicate mushy zone. Dimensions of domain are in m. 

Fig. 6. Simulated results of Al-2 wt% Cu at 2500 s with lateral temperature gradient of the 
domain width of 6 mm: (a) temperature field; (b) thermosolutal convection and contour lines 
indicating mushy zone. Dimensions of domain are in m. 
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Fig. 1. Typical morphology of liquid-solid interface during plane front directional solidification 
of (a) Pb-2.2 wt% Sb and (b) Al-2.2 wt% Cu alloy. Notice the convex liquid-solid interface for 
the Al-Cu alloy as compared with nearly flat interface for the Pb-Sb alloy. The radial solute 
segregation has lead to the break-down of the planar front near the periphery of the Al-Cu 
sample, where cellular morphology has evolved. 
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Fig. 2. Macrosegregation along directionally solidified length of Pb-2.2 wt% Sb alloy grown in 
7 mm ID quartz crucibles at 0.4 pm s' 1 and 120 K cm' 1 ; and of an Al-2 wt% Cu alloy grown in 6 
mm ID alumina crucibles at 90 K cm' 1 and 0.5 jitm s' 1 . Each alloy was grown in the presence of a 
transverse magnetic field of 0.4 T (MF) and without it (NMF). 
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Fig. 3. Influence of reducing the ampoule diameter on the longitudinal macrosegregation in 
directionally solidified Pb- 2.2 wt% Sb alloy samples (plane front). The 1 , 2, and 3 mm 
samples were grown at 0.4 pm s' 1 , 82 K cm' 1 , while the 7 mm sample was grown at 0.4 pm s' 1 , 
120 K cm' 1 . 
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Fig. 6. Simulated results of Al-2 wt% Cu at 2500 s with lateral temperature gradient of the 
domain width of 6 mm: (a) temperature field; (b) thermosolutal convection and contour lines 
indicating mushy zone. Dimensions of domain are in m. 



